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ABSTRACT The importance of transition modeling in the computation of 
compressible, unsteady separated flows is discussed. The study showed that 
it is critical to predict the experimentally attained transition point properly 
in order to obtain good agreement with data at the same Mach number and 
Reynolds number. 


1. Introduction 

Recent advances in computing power have made it possible to obtain solutions of 
the 2- and 3-dimensional Navier-Stokes equations for complex separated flow problems, 
such as dynamic stall on airfoils and finite-span wings. Calculations by Clarkson et. 
al[l] have shown that simple algebraic and half-equation turbulence models are inade- 
quate to predict the measured dynamic stall hysteresis loops. It is well recognized that 
leading-edge stall is preceded by the development of a. leading-edge separation bubble 
for chord-based Reynolds number < 1 x 10 6 . The size of the bubble is determined by the 
transition and entrainment processes within the bubble. Two key parameters, namely 
the transition onset location and the transition length need to be modeled correctly for 
successful prediction of the bubble and its bursting, as shown in Refs. [2]- [4]. Com- 
pressible dynamic stall studies by Ekaterinaris and Platzer[5] at Re = 4 x 10 6 showed 
that transition plays a dominant role even at full-scale Reynolds numbers since it is a 
leading-edge type of stall and hence, the physics of transition needs to be incorporated 
in modeling the flow. 

In the present paper, the deep dynamic stall flow over a NACA 0012 airfoil is studied 
at Re = 1.1 x 10 6 , where corresponding experimental data are available[6]. However, the 
transition modeling used in Ref. 5 for high Reynolds numbers is replaced by that used by 
Ekaterinaris et al[2] for transitional Reynolds number flows. The methodology is based 
on use of Michel’s criterion for determination of the transition onset location, (given 
by the matching of the two functions Reg(x) and f(x) = 1.174(1.0 + 22400)/i?e°- 16 ), 
with Re z and Reg based on u 0 0 rather than the boundary layer edge velocity u e . The 
transition length is obtained using the Chen-Thyson model: 
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with G- /tr = 213[/o</(7?e Jtr ) - 4.7323]/3 as the transition constant. 
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2. Experimental Investigations of Transition Effects on Dynamic Stall 

In [6], the Reynolds number effects on oscillating airfoil compressible dynamic stall 
flow have been investigated on two NACA 0012 airfoil models at constant Mach Num- 
bers, with and without tripping the boundary layer on both models. The corresponding 
(untripped) values of Rc r were 0.54 x 10 6 and 1.1 x 10 6 respectively. The experiments 
showed that with increasing Reynolds number, the length of the laminar separation 
bubble was reduced. Further, dynamic stall onset was shifted to higher incidences dur- 
ing the upstroke of the airfoil and the flow was able to withstand higher adverse pressure 
gradients. In the present study different transition onset models are applied to simulate 
the experimental conditions, i.e. natural transition and boundary layer tripping. 

3. Implementation of Transition Model into Navier Stokes Code 

The present numerical study uses a time-accurate (Beam and Warming) 2D Navier- 
Stokes code[7]. Calculations until now used fully turbulent flow during the entire os- 
cillation cycle. As shown in [6], since transition plays a dominant role in the initiation 
and further development of the dynamic stall process, it has to be taken into account 
in numerical modeling as well. It is even more important to account for transition if 
research towards favorably influencing dynamic stall (or other aerodynamic properties) 
is pursued by means of dynamic airfoil shape modifications. Recent numerical studies 
by Geissler and Raffcl[8] and Geissler and Sobieczky[9] as well as wind tunnel measure- 
ments by Chandrasekhara et al [10] have already shown the considerable benefits of such 
an approach. 

In the present numerical study, different measures to determine transition onset 
have been investigated on a NACA 0012 airfoil section in oscillatory pitching motion. 
To match experimental data in [6] the following set of parameters have been used in the 
study: incidence, a = 10° + 10°sinu!<; Mach number, M = 0.3; Reynolds number, Re = 
1.1 x 10 6 ; reduced frequency , k = = 0.05. For the determination of transition onset, 

three different options have been used: (a) prescribed transition-onset corresponding to 
the tripping device in the experiment; (b) calculation of transition onset by Michels 
criterion; (c) calculation of transition onset at the instantaneous position of the pressure 
minimum. 

4. Results and Discussion 

Figure 1 shows a typical result of the present investigation. The pressure minimum 
is plotted versus incidence during the airfoil upstroke. The solid curve is obtained 
if fully turbulent flow is assumed. Dynamic stall onset is characterized by the rapid 
development of the dynamic stall vortex which occurs approximately at the maximum 
of this curve. 

The dotted curve is obtained by prescribing the transition onset at x/c = 0.015 to 
match the leading edge of the transition strip (located between x/c = 0.015 and 0.030) in 
the experiments[Gj. Remarkable agreement of this (dotted) curve with the experimental 
(tripped) data can be observed. Changing transition onset upstream of x/c = 0.015 to 
x/c=0.005, (not shown) shifts the maximum of the curve towards the fully turbulent 
limit. Shifting transition onset downstream towards the end of the transition strip at 
x/c = 0.03 leads to an earlier dynamic stall onset (at a = 12 deg) and shows a trend 
towards numerical instability. 

The dashed curve in Fig. 1 indicates the result obtained with Michel's criterion 
(assumed to approximately determine natural transition onset), which has been proven 
to be applicable for Re > 1 x 10 6 . Also shown in the plot is the experimental result 
for free (untripped) transition. The maximum of the experimental curve is obtained 
prior to the maximum obtained with Michel’s criterion. In both experimental and 
numerical cases a trend to shift dynamic stall onset slightly beyond the point obtained 



with tripping can be observed. 

The long-dashed curve in Fig. 1 finally indicates the case of transition onset fixed 
to the instantaneous position of the pressure minimum, which is located very close to 
the airfoil leading edge. This results in a delay of dynamic stall onset towards the fully 
turbulent case. 

Figure 2 compares the calculated positions of transition onset as obtained by 
Michel’s criterion (dashed curve) and from the C» min locations (solid curve). Also 
shown are the fixed positions of transition onset wliich have been used in the present 
calculations. The pressure minimum moves rapidly upstream from about 14% chord at 
a — 0 deg. to xjc « 0.005 for a — 10 deg. Michel’s criterion predicts transition onset 
further downstream and reaches the 0.5% location beyond a = 15 deg. Indicated in the 
plot are also the (constant) positions of fixed transition onset. The distance between 
these lines and the pressure minimum indicates the area where the laminar boundary 
layer is subjected to an adverse pressure gradient. 

Figure 3 shows the calculated force and moment coefficient loops for the cycle. 
Two different cases are illustrated: the fully turbulent flow (solid curves - the reference 
case) and the case of transition onset fixed at x/c = 0.015. A considerable reduction 
of maximum lift can be observed in the transitional case. Dynamic stall onset occurs 
about 4-5 deg earlier within the oscillatory loop when transition modeling is included. 
Similarly, the drag-rise and the beginning of moment stall occur earlier in the cycle. 

5. Conclusions and Future Investigations 

With the implementation of a transition model into the time-accurate Navier-Stokes 
code the critical role of transition on dynamic stall onset has been documented. The 
rather simple Chen-Thyson model which has been used in a quasi-steady manner gives 
results comparing reasonably with the experimental data. The study has also addressed 
the problem of determination of the transition onset location. If tripping is used in 
the experiment, best agreement with numerical calculations is obtained by fixing the 
transition onset point at leading edge of the transition strip. Michel’s criterion gives 
transition onset locations which seem to be too far downstream. Nevertheless, the 
results of the present study show trends similar to those seen in experiments. Use of 
suction peak location for transition onset pushes dynamic stall onset towards the fully 
turbulent limit. The present results provide useful design guidance for the modification 
of the airfoil shape to delay dynamic stall onset. 
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